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The formation of chain-folded structures from the melt is 
observed in molecular dynamics simulations resembling the 
lamellae of polymer crystals. Crystallization and subsequent 
melting temperatures are related linearly to the inverse lamel- 
lar thickness. Analysis of the single chain conformations in 
the crystal shows that most chains reenter the same lamella by 
tight backfolds. Simulations are performed with a mesoscopic 
bead-spring model including a specific angle bending poten- 
tial. They demonstrate that chain stiffness alone, without an 
attractive inter-particle potential, is a sufficient driving force 
for the formation of chain-folded lamellae. 



For polymeric materials, the connectivity of the long- 
chain molecules usually hinders the formation of perfect 
crystals. Instead, most crystallization conditions lead to 
partially crystalline structures consisting of a stacking 
of chain- folded lamellae and amorphous regions. There 
is nowadays wide agreement on the structure of this 
semicrystalline phase jj]. However, the pathway which 
leads to it and why it is taken is one of the unsettled ques- 
tions in polymer science. Actively debated ideas (^-0] 
are intermediate mesomorphic phases g] and the possi- 
bility of a spinodal mode in contrast to nucleation 
and growth. Existing theories of polymer crystalliza- 
tion jH)|-|l2| are still considered incomplete 
Thus, it is desirable to develop new tools for examining 
nucleation and crystal growth, and especially its early 
stage. 

Computer simulations generally have two advantages 
which may help to understand the noncquilibrium pro- 
cess of polymer crystallization. First, the trajectories 
of all particles are available. This can give information 
about single molecules which are not accessible by ex- 
periments. Second, models of variable complexity can 
be simulated. The use of successive levels of approxi- 
mations has the advantage to give insight into the ef- 
fect of certain interactions separated from others. In 
this spirit, this Letter explains how chain- folded struc- 
tures emerge from the supercooled melt when only the 
excluded volume and the possible conformations of the 
polymer chains are taken into account. 

Many other simulation approaches have already been 
proposed in the context of polymer crystallization. 
Monte Carlo (MC) techniques were used to test theories 
of the lamellar thickness with lattice models |Ii|,[l^.|ig| . 
On a larger length scale, MC was successfully applied 



to model 2D crystallization |17| ]. Molecular dynamics 
(MD) is often more realistic, but simulations have been 
restricted either to few chains and short times [18], or 
to structure formation of very short chains Jl9|— 21 1 , or 
to an isolated long chain which collapses into a "crys- 
talline" globule |2^]. Recent MD simulation studies 
which come closest to our work deal with the crystal- 
lization of 20mers from the melt |fl9| , of lOmers in thin 
films [^TJ, and of few long chains from solution p3[ . Here, 
we present large scale MD simulations of structure forma- 
tion from the melt for a wide range of chain lengths and 
temperatures. We examine the homogeneous nucleation 
regime, since no walls or starting nuclei are introduced. 
In this Letter, we demonstrate in particular that systems 
of chains long enough to fold can be simulated from the 
amorphous melt up to the chain-folded lamella. The in- 
dividual particle trajectories are then used to analyze the 
nucleation process and the chain folding in detail. 

We use a simplified version of a model for poly(vinyl 
alcohol) (PVA) derived by a systematic coarse-graining 
procedure p4| from fully atomistic simulations. A coarse- 
grained bead represents one monomeric unit of PVA. (In 
contrast, the often used united-atom models 19, 20|,|23|] 
consider one bead per carbon.) The beads are connected 
by harmonic springs, and additionally interact by an an- 
gle bending potential retaining information on the tor- 
sional states of the atomistic backbone. The nonbonded 
interactions are approximated by a Lennard-Jones 6-9 
potential. Furthermore, for the simulation of a dense 
melt, it is sufficient to take only the repulsive part of the 
potential f25|] . This leads to some shift of energy and 
pressure while it affects only little the qualitative behav- 
ior of the system. This approximation speeds up the 
simulation. But it also switches off one possible driving 
force of crystallization, the cohesive energy, which would 
increase with decreasing temperature. Thus, when low- 
ering temperature, only the chain stiffness will increase. 

Our model is defined by several parameters p6[ . 
Length scales are fixed by the mapping from atomistic 
simulations. Units are reported in a = 0.52 nm cor- 
responding roughly to the chain diameter of PVA. The 
bond length is 0.5a = 0.26 nm. In these units, the density 
is about 2.1 monomers per a 3 . The nonbonded poten- 
tial between monomers is V(r)/ks = — 0.75/r 6 + 0.53/r 9 
which is cut and shifted to zero at the minimum r cut = 
1.02cr. Nonbonded interactions between first and second 
neighbors along each chain are excluded. The angle po- 
tential shown in Fig. |l| has three minima at 180°, 126°, 
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FIG. 1. Sketch of the coarse-grained model "CG-PVA". 
Each bead represents a monomer with two carbon atoms on 
the backbone. Successive beads are connected by harmonic 
springs. Three successive beads interact via the shown angle 
bending potential. 

and 95° corresponding to trans-trans, trans-gauche, and 
gauche-gauche conformations of the atomistic backbone 
chain. Note that this model, though initially fitted to 
represent PVA, is also suitable for polyethylene which 
has almost the same dimensions as PVA. The equations 
of motion are integrated by the velocity- Verlet algorithm 
using a timestep of 0.005 r |27| . Temperature is kept con- 
stant through a Langevin-thermostat with friction con- 
stant T = 0.5. Temperatures and energies are expressed 
in dimensionless units with m = ks = 1; T = 1 corre- 
sponds to the high temperature phase of the amorphous 
melt. Constant pressure is imposed by the isotropic 
Berendsen manostat |2§| l . The simulations are performed 
in a cubic box with periodic boundary conditions. 

The model is studied for chain lengths N = 10, 20, 
50, and 100 by two kinds of experiments. The first con- 
sists in cooling at constant rate, the second is isothermal 
relaxation: the system is suddenly quenched to a temper- 
ature T c below the melting temperature. In both cases 
we start with an amorphous configuration equilibrated at 
T = 1. Figure ^| shows the volume as a function of tem- 
perature in cooling and heating cycles. It decreases with 
temperature; a sudden drop of the volume is attributed 
to a phase transformation. Note the hysteresis. Indeed, 
crystallization only takes place when a stable nucleus is 
formed, i.e. at a certain supercooling below the melting 
temperature. In addition, the hysteresis is enlarged due 
to the fast cooling and heating rates. Once formed, a 
crystal is also stabilized by the periodic boundary condi- 
tions which lead to a higher melting temperature. 

Figure ^| also shows that shorter chains begin to crys- 
tallize at lower temperatures than longer chains. Simi- 
larly, the decamers melt at lower T m than the icosamers. 
However, still longer chains melt at about the same tem- 
perature as the icosamers. Furthermore, the volume 
change for N = 50 and N — 100 is less pronounced 
than for N — 20. Both observations indicate that the 
crystallinity of these samples is far from 100%. In fact, 
the shorter chains N = 10 and N = 20 always form ex- 
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FIG. 2. Volume dependence in cooling/heating cycles for 
different chain lengths. The cooling rate is 5 x 10 _6 r _1 , the 
heating rate 2 x 10 _5 r _1 . 

tended chain crystals. At temperatures below T w 0.7, 
crystallization sets in immediately and small domains are 
formed similar to those reported in Ref. |]l9| . For smaller 
temperature jumps, the formation of a single nucleus is 
observed which grows in regular layers. More details will 
be given elsewhere . 

Let us focus on the longest chains, N = 100, in the 
following. They always form folded lamellar structures 
(see Fig. ||). The stems (the completely stretched parts) 
are much longer at higher temperature. Furthermore, the 
structure factor (not shown) indicates long range order 
with Bragg peaks of the hexagonal lattice. The structure 
factor also has a peak at small wave vectors correspond- 
ing to the thickness of the crystal lamella. In Ref. p(| 
we discuss an alternative way of determining the lamellar 
thickness with the bond orientation correlation function. 
It gives an upper bound of the stem length. Since most 
chains form tight backfolds (see below), this is a good 
estimate of the lamellar thickness. Figure ^ shows that 
the inverse of this average stem length increases linearly 
with decreasing T c , a relation observed experimentally 
for most crystallizing polymers. The same holds for the 
melting temperature of the corresponding structures. At 
T c = 0.7, most chains form 3 folds, at T c = 0.74 they have 
mostly 2 folds. At T c = 0.78 the first nucleus grows with 
once-folded chains; however, it continues to grow with 
a mixture of once and twice folded chains which do not 
unfold during the simulation time. At T c = 0.8, finally, 
most chains exhibit only one fold. 

This can be summarized as follows: the higher the tem- 
perature the better the crystal. Though perhaps para- 
doxical at first sight, this is generally observed for crys- 
tallizing systems and in particular for polymers because 
the thermodynamically favorable extended chain crystal 
is kinetically not accessible. For a single chain, crystal- 
lizing with many folds is much closer to the initial coil 
state in the melt than a once folded or even extended 
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FIG. 3. (color) Two examples of structures obtained at 
temperatures T c = 0.7 (a) and T c = 0.8 (b) with 192 chains 
of TV = 100 monomers. Note the larger lamellar thickness 
at T c = 0.8. At T c — 0.7, two smaller domains with differ- 
ent orientation are formed. Color coding is according to the 
bond angle (cf. Fig. |l|): blue=tt, red=gt, yellow=gg. The 
box length is about 10 nm. Note that periodic boundary con- 
ditions are applied; chains are not cut but mapped with their 
center of mass into the central box. 

chain. So, at lower temperatures, when the mobility is 
lower and the nucleation probability is higher, the chains 
will end up in crystals of short stems before having time 
to unfold. The barrier between crystals of different fold 
lengths is high and is overcome only with additional ac- 
tivation, for example, by annealing. 

All simulation results reported so far are consistent 
with the interpretation of well known experimental find- 
ings. This gives us confidence in the model proposed. 
On this basis, we may use the simulation as a nanoscope 
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FIG. 4. Dependence of lamellar thickness on the crystal- 
lization temperature and melting temperature of the corre- 
sponding final structures. The lines are linear fits. Melting 
temperatures are rather overestimated due to finite size ef- 
fects. 
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FIG. 5. (color) Snapshots of the crystallization at T c — 0.8. 
(a) The 13 chains participating in the initial nucleus are shown 
(blue). The following snapshots contain 20 more chains (thin 
lines) which align to the crystal in the first 40000r. (b) The 
red chain is partly attaching to the crystal front, (c) The 
hairpin of this chain is growing to fill the whole lamellar thick- 
ness, (d) Finally, the hairpin is vanishing: the chain glides 
until having only two stems in the crystal and a large bridge 
on the fold surface. 

and visualize the evolution of single chains as well as the 
conformations of the chains in the crystal. Figure ^ il- 
lustrates snapshots from the crystallization at T c = 0.8. 
In this case a slightly different temperature protocol was 
applied: the start nucleus was obtained by quenching 
the melt from T = 1.0 to T = 0.7 over 10000 t. Then, 
the temperature was raised to T = 0.8 which ultimately 
leads to the structure in Fig. ||b after 140000 r. Initially, 
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the nucleus grows primarily in the direction of the stems 
(a,b). Later on, further stems reorient and "attach" to 
the primary nucleus. This "attachment" can start with 
arbitrary parts of the chain. However, the chains move 
until aligning their ends with the existing crystal. In the 
sequence of Fig. || one chain is highlighted which joins 
the crystal in a complicated way. It attaches first via a 
hairpin which advances until aligning with the total stem 
length of the crystal (c). However, the other end of this 
chain is quite far. The chain slides along its contour to 
finally annihilate the hairpin. 

A quantitative analysis of the single chain conforma- 
tions results in 3/4 of the folds in the final structure at 
T c = 0.8 being backfolds. One third of them are hairpin 
folds, only 1/4 of the backfolds reenter the crystal farther 
than at the 3rd neighbor position. At lower T c , the frac- 
tion of backfolds in our samples was found to be around 
60%. Due to domain formation (Fig. ||a), up to 25% of 
the folds link stems participating in different domains. 
The distribution of the re-entry distances of backfolds is 
almost the same for all crystallization temperatures. 

In summary, this Letter highlights results of a versa- 
tile model suitable for studying the chain-folding process 
of polymer crystallization. The chosen level of coarse- 
graining retains enough information on the local confor- 
mations while the speed-up due to the reduction of de- 
tails is sufficient to follow nucleation and growth from 
the melt up to chain-folded structures. Starting from 
the chemistry of PVA, the model finally consists only of 
excluded volume interaction and a specific angle bend- 
ing potential favoring stretched conformations and al- 
lowing precise fold states; there is no explicit attraction 
between monomers. We understand the crystallization 
in our model as driven by an energetic and an entropic 
component. When lowering temperature, energy leads 
through the angle potential (Fig. [[]) to an increased chain 
stiffness (the persistence length of an isolated, free chain 
passes from 3.4 monomers at T = 1.0 to 6 monomers 
at T = 0.7). Entropy favors the parallel ordering as 
in liquid crystal j29], because the gain in translational 
entropy of the ordered structure at lower temperature 
is possibly much larger than the loss in conformational 
entropy. The temperature dependent kinetics then de- 
termines the form of the nonequilibrium structures of 
chain-folded lamellae. Yet, the qualitative chain length 
dependence of crystallization and melting temperatures 
(Fig. ||) as well as the lamellar thickness behavior (Fig. ||) 
match well experimental findings. This underlines that 
the orientational and conformational ordering play a ma- 
jor role for polymer crystallization. It further suggests 
that crystallization might pass through an intermediate 
step of liquid crystalline order before locking in into the 
final crystal structure [jfj, at least in a region in front of 
the growth front. The model could now be used to ad- 
dress questions of recent crystallization experiments in 
confined geometries and, with larger systems, the open 



question of a spinodal initiation of crystallization. 
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